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Recently Maleknejad and Sheikh-Jabbari have proposed a new model for inflation with non- 
Abelian gauge fields (Gauge-flation), and they have studied the model by numerical methods Q|. In 
this model, the isotropy of space-time is recovered by suitable combination of gauge configurations, 
and a scalar field is constructed by gauge field and the scale factor, which produces inflation period. 
In this work, exact solutions for the scalar field and the Hubble parameter are presented and we 
provide analytic solutions for the numerical results. We explicitly present Hubble parameter and 
fields as functions of time and it is also demonstrated that in some conditions they are damped 
oscillator. Moreover, reheating period in the model is discussed. 



(N 

o 

Or 
OS 



cr 

i ' 

Sh ' 



in 
> 
o 
in 

o 

(N 



X 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Observations of the cosmic microwave background and 
large scale structure are consistent with the isotropic 
symmetry of space-time and inflation a period of 
accelerated expansion occurring in early universe. To 
save isotropic symmetry of background space-time, many 
models use a single or multi-scalar field(s)(inflaton) with 
ansatz that the inflaton field rolls slowly down its po- 
tential (inflationary period) Eventually the field(s) 
oscillates around the minimum of its potential and de- 
cays into light particles (reheating period). It has been 
assumed that reheating period took place just after in- 
flationary period and the field behaved like dust matter 
in the reheating period [l[. 

We have had few models that can be solved exactly in 
the cosmological context. The well-known example, that 
can be solved exactly, is the exponential potential P, 0] 
which, with suitable parameters, gives us inflationary so- 
lution(power law inflation). Usually, solving field equa- 
tions are limited to either the inflationary period, with 
slow roll approximations, or reheating period. 
The success of gauge symmetry in the standard parti- 
cle physics shows that the gauge symmetry is the cor- 
rect tool to construct an effective field theory. But, at 
first glance, it seems impossible to reconcile the isotropic 
symmetry with the vectorial nature of gauge fields. How- 
ever, recently a new model for inflation, Gauge-flation, 
based on gauge field , A a ^, where a, b, . . . and fi,v,... 
are used for the indices of gauge algebra and the space- 
time respectively, has been proposed by Maleknejad and 
Sheikh- Jabbari(MS) 0. The rotational symmetry in 3d 
space is retained by introducing three gauge fields, such 
that these gauge fields rotate among each other by SU(2) 
non-Abelian gauge transformations Q. In this model, a 
scalar field that causes inflation period, ip, is constructed 
from the gauge field and the scale factor (see below). MS 
have studied the model by numerical analysis and ob- 
tained numerical solutions for the model Q. 
In this work we investigate the model by analytic method. 
As we will see, the model has some features that allow us 
to obtain leading order of fields in all epochs of the uni- 



verse and not just only in a specific period. We provide 
an analytic expression for MS solutions and we study re- 
heating period in this model. 

This letter is organized as follows: in §11 we briefly re- 
view the model and obtain equations in terms of variables 
which we will use in this letter. In §111 we consider a spe- 
cial solution in spectrum of solutions of the model, that 
the equations can be solved exactly. The solution yields 
analytic expression for the Hubble parameter and fields 
for all epochs of the universe. Also we give the general 
form of solutions for the fields. In §IV we obtain the 
analytic form for the MS solutions. Constraints on the 
solutions from slow roll conditions and observations are 
discussed in §V. In § VI we discuss about reheating pe- 
riod in the model. We summarize our finding in §VII. In 
Appendix A, we obtain the equation of motion of fields 
in terms of variables that we use in this letter. 



II. THE MODEL 

We work with a general flat-space FRW background 
metric with signature (— -f_H — h), and reduced Plank 
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Lagrangian, 

L = 
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Following Q, we consider the effective 
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where e^ v is the totally antisymmetric tensor and the 
strength field F is 
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where e a b c is the totally antisymmetric tensor. Sheikh- 
Jabbari [if has shown how F 4 term in ([lj is obtained by 
integrating out a massive axion field in Chromo— Natural 
inflation [5(. Also, it has been argued how other dimen- 
sion 8 level terms are suppressed [J] . 
To retain isotropy symmetry of space-time we have to set 
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Applying Einstein's equations with this setup results in 
I 



H + 2H 2 = k 2 - 



(4) 
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where a is the scale factor. The equation of motion for 
4> is obtained by combination of the equations in((3]) Q, 
so it is not an independent equation and we take (|4]) 
to study the model . Neither 4>{t) nor a{t) are scalar 
under general coordinate transformation, while ip(t) = 
4>{t)/a{t) is a scalar, so we seek solutions for it. We define 
a new variable a(t) by writing 



— = vMO cosa(i), 
^Tjj = vW) sina(i), 



(5) 



where e(t) = —H/H 2 and a(t) is an arbitrary function of 
time. Using ([5]), the system of equations in (QJ is reduced 
to 



e(t) 



-e 2 {t)H 2 sin 2 2a(t), 



(6) 



noting that < e(t) < 2. Without assumption either 
on the H or on the a(t) Eq. ([5]) cannot be solved to 
obtain an expression for the Hubble parameter as a func- 
tion of time(except for k = that results in H = l/2t, 
" radiation epoch" ) . To explore the model, we define 



Hsm2a(t) = fif(t), 



(7) 



where f3 is an arbitrary number (dimcnsionless in Plank 
units), but P < H. 

As we will see, the MS solutions can be obtained by a 
special form for the f(t). Before we study the MS so- 
lution we take a "simple" ansatz for f(t), that not only 
helps us to understand the MS solutions, but also, it gives 
information about dynamics of the model. 



III. THE SIMPLE ANSATZ 

In this section we use the following simple form for f(t) 

f(t)=smut, (8) 

where ui is an arbitrary number. Using Eqs. (j6|) and ([8]). 
we have 



\/l + 2fcsin 2 uxt- 1 



| sin 2 uit 



(9) 



where k = (k(3) . Note that lira e(t) = 2, where n is 

ujt— >mr 

an integer number. e(t) is a periodic function of time and 



for k > 4 it crosses e(t) = 1 line (acceleration expansion 
phase), and for k > 8 x 10 4 it crosses e(t) = 10~ 2 line. 
The e(t) is designed as in Fig. 1. 

To understand behaviour of e(t), we can derive an useful 
formula for e(t) in inflationary period(i.e., when e(t) -C 
1). From Eqs. j6]) and ([7|), e(t) in this limit becomes 



e(t) 



V8 



Vk sin tot 



(10) 



inf 



where the subscript "inf" denotes that the above 
expression is valid when e(t) -C 1, i.e., < wti n f < it. 
For large value for k, the validity of Eq. (fTU)) is broken 
when t is very close to mr/u), so the universe almost 
evolves through acceleration expansion phases for fc ^ 1 , 
but eventually exits abruptly from inflationary period, 
as indicated in Fig. 1. 




FIG. 1. The e(t) versus t for k = 10°, u = 2. After inflation, 
the universe exits abruptly from inflationary period. 
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FIG. 2. The Hubble parameter versus t for k = 10 , to = 
2. Since —H increases with time, the second steep slope is 
shorter than the first steep slope. 
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FIG. 3. The Hubble parameter versus t for k = 10 , ui = 2. 
Integration of ((9]) gives 

H(t) = ku> 1 2 ( {\ - y/l + k(l - cos2wi)) cot cut + 



{l + 2k)F[ujt\-2k] -E[cot\-2k] 



(11) 



where F[iot \ — 2k] and E[u>t \ — 2k] are elliptic integrals 
of the first and second kind respectively and the specific 
combination of them in (|lip is increased monotonically 
with time 01. Also 



lim ( (— 



((-1 + \j\ + fc(l - cos2wi)) cotwt = 0, (12) 

so the Hubble parameter is damped (but it is singular at 
t = 0) . The Hubble parameter is planned as in Fig. 2 and 
Fig. 3, where it has gentle slopes when e < 1( inflationary 
phase). Since H is decreased with time and e(t) is a 
periodic function of time, we conclude that not only H < 
but also — H is decreased with time, as indicated in Fig. 
2. 

If we demand just one inflationary phase for the universe 
and define At in f as the duration of time that inflation 
takes place, from (fT0|) we have 



jAtinf < 7T. 



(13) 



Note that (|13|) is the upper bound on uiAti n f that the 
model predicts itself and it must be consistent with ob- 
servational constraints (see below). 
From Eqs.([5]) and (J7J), we find 



1± a/1 



P 2 sin- cut 



H 2 



fJ 
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where e(t) and H are given by ((9]) and (fTTj) respectively. 
It is worth to mention that our results in (j9|), (fTTj) and 
(fT4"| arc valid for all values of e(t) and this is one of 



the interesting properties of the model compared with 
other models for inflation, which use e(t) as a perturba- 
tion parameter and their solutions for the Hubble rate or 
field(s) are limited to a specific epoch (inflationary period 
or reheating period). In this sense, the above expressions 
are exact (non-perturbative) and nonsingular solutions of 
nonlinear equations (|4|). 

The inverse proportionality of ip± to ^fg, shows that they 
are not exist in perturbation regime. 
To understand the behaviour of (fl"4"| , we use the fact 
that j3 < H, and expand (|T4"]) that gives 



fJ 



2<l 



sin ut. 



(15) 

Here the higher terms ~ (3 2 / H 2 have been neglected. So, 
ip+ is damped oscillator and ip- is oscillator (till j3 < H). 
The leading order behaviour of ip± are sketched in Fig. 
4. 

Recall that in the reduced Plank units, H < 1( after 
Plank time), and /3 < H, so ip- is smaller than tp+ in the 
inflation period. However, from the Friedman equation, 



3H , we realise that both of them have the same 



density, and the behaviour of H (in Eq. (|TT|) and Fig. 
3) shows that the cosmic expansion dilutes the density 
of fields, although the expectation value of tp+ is greater 
than ip- after inflation. 

For general type of f(t) in (JT)), H cannot be obtained in 
terms of well-known functions, but by algebraic manip- 
ulations of equations JS]),© and (O, one can show that 
for leading order of fields we have 



'/(*). 
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(16) 

where Hf is the Hubble rate for f(t). Since Hf is usually 
unknown, the above expression is formal for ip + f , but 
the leading order ofip—f is given by (|16[) for any f(t). 
Note that Eqs. (|T5)) are valid for all values of e(t). The 
physical meaning of fit) is clear from (p~6|) . i.e., fit) 
shows boundary conditions on the fields. 



IV. THE MS SOLUTIONS 

We pointed out that MS have analysed the model by 
numerical methods. Here we will show that to rederive 
the results we must take the following form for f(t) 

m = ( M<*t + cf + ---) \ 2 sfMg(Jt)t (17) 



ai t 



e(t) and the leading order of ip- are obtained by ([5]), ©, 
(IT61) and JIT]) as 



2 (yi + 2kfi Is it) - 1 



(18) 
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FIG. 4. The leading order behaviour of ip_ (solid line) and 
ip\ (dashed line) for k = 10 6 , u = 2, /3 = g = 1(T 6 



Figures 5 and 6 are obtained by our analytic method and 
they have the same pattern as MS have obtained^. 
Note that if /ms —> then e(t) = 2. Although in this case 
we cannot obtain an exact expression for the Hubble pa- 
rameter (for all time), but with (fT7|) . we can obtain some 
information about the MS solution in various regimes. 
The most important coefficient of (|I7|) in inflationary 
phase is a.\. When ot\t < 1 we have f(t) w 1, so in 
this regime from ((T7)) and (US]), we have 



£Ms{t) 



2(x/l + 2fc - 1) 
k 



(19) 



i.e., 6Ms(t) is constant in this regime as indicated in Fig. 
5. To have e(t) < ICT 2 , we must take k > 8 x 10 4 . 
Therefore the Hubble parameter in this regime is 



H M s(t) 



2(VT + 2k- l)t 
Hence, from (|T6| and (fT9|) it follows that 



(20) 



The Hubble parameter in this regime can be derived by 
integration of (f2"2"|) . that is 



1 



(23) 



similarly, the leading order of ip± are given by the fol- 
lowing formulas 



i'-Ms(t) 



2gt 
2g7r 2 



a\k 

(t--r 

ax 



(24) 



for t > 1, the other terms, (a2, ■ ■ • ), in (|I7|) arc impor- 
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FIG. 5. The e(<) versus t for = 10 6 , cj = 2, ai = 10,q 2 
0,«3 = 2. Compare with the figures in 0] 



tant, and fMs(t) oscillates faster, so e(t) and ip± oscillate 
faster. 

Finally, for t 3> 1, the denominator in (|I7p is increased, 
and results in /j\/s ~ 0. Therefore in this regime, 
e(t) — > 2, ip_ — > as indicated in Fig. 5 and Fig. 6. 



2 5 2 ( v / TT2fc- 1) *' 



(21) 



^_ / 2(yTT2fc- 1) 
2.9V * 



so, the leading order of ip— is constant in this regime, as 
indicated in Fig. 6. 

When ait 7r, the numerator in ()17|) equals to zero, and 
therefore cms ~> 2. Just like the simple ansatz in §11, 
for large k, the universe almost evolves through inflation 
period but eventually exits abruptly when t — > n/ai. 
Using (|17p . the Taylor expansion of e(t) in p^|) around 
t = n/ai is 



k 



(t 



ai 



(22) 
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FIG. 6. The ip- versus t with the same parameters as in Fig. 
1 and P = g = 10 -6 . Compare with the figures in 0] 
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Also, using (jT6|) . wc have 

H MS {t) -> — , V- 



1 

2f 



1 1 



V2g t 



(25) 



Results from various numerical analysis of the model 
with various initial values Q have similar properties of 
the MS solutions, that lead us to the following conjecture 
about the model: the MS solution, (|17p . is the attractor 
solution of the model. 



V. COMPLEMENTARY SLOW ROLL 
CONDITIONS 

If we demand that our fields be inflaton fields, comple- 
mentary slow roll conditions are required in inflationary 
period ,i.e., not only e(t) <C 1 but also 



H f S. 



< 1. 



f°±f 



(26) 



For ip+fi the conditions do not have any restriction on 
the parameters of rf+f, therefore we will focus on ip-f- 



A. Complementary slow roll conditions for the 
simple ansatz 

For the simple ansatz, ([8]), the conditions in (|26| yield 

uj sin 2cuti n f 



H 



<C 2 tanwtm/ 



H« 2 



(27) 



If < cutinf < 7r/2, the conditions in ([27]) are reduced to 
ui <C H , that is agreement with (flU|) . 
If we take uiAti n f = 1/2, then using the condition uj <C 
H , wc obtain 



HAt mf > 1. 



(28) 



The current cosmic microwave background data indicate 
that during inflation epoch, H < lO~ 5 M pl Q. If H m 
10~ 5 M p i, the relation (p8|). implies that At inf > 10%;. 
Number of e-folding is given by 



tf 

-- J Halt, 



(29) 



If we use uj Atinf = 1/2, then assuming that tf > 10ij, 
the numerical integration of the Hubble parameter 
in (jlip shows that to have N e > 60, we must take 
k > 4 x 10 6 . 

So, if we set /3 = 10~ 6 ( in Plank units), then to have 
sufficient e-folds, we must take k>2x 10 9 . 



B. Complementary slow roll conditions for the MS 
solution 



As for the MS solutions, ([5BJ yields 



2(aiicot(aii) - 1) < H M s t , 
1 — (ait cscctit) 2 



T + OL\t cot oc\t 



•C Hms t 



(30) 



Here we have neglected other terms that do not have any 
effect in inflation period. For ait < 1 we have 



2(ai*cot(aii) - 1) 
1 — (ait cscaii) 2 



2(aitf 
3 
2 



(31) 



1 + T^(ait) 2 ). 



Using Eqs. (|20l) . (|3U|) and (|31l) . the complementary slow 
roll conditions arc reduced to 



2<yi + 2k - 1; 



»1, 



(32) 



hence, it is sufficient that fc > 8 x 10 4 . 
Just like the simple ansatz, if we take aiAti n f = 1/2, 
from relations (|30[) and (f3~Tj) . we have HAU n f 3> 1. 
Number of e-folding is given by (|20]) and (|29|l as 



In 



2(Vl + 2fc- 1) W 
so, for fc > 8 x 10 4 we have 

tf 



N e > 100 In ■ 



(33) 



(34) 



The general cosmological perturbation of the model was 
developed in Q. 



VI. 



REHEATING 



In the most models for inflation to have successful re- 
heating period, it is necessary that after inflation period, 
inflaton(s) behaves like dust matters, and then decays 
into relativistic matter. But in the Gauge-fiation model 
the fields can be decayed into relativistic matter, with- 
out going to dust matter phase as we saw in §IV. One 
way to see this point is to see the Lagrangian ([1]), in the 
inflationary period the (FF) 2 term is dominate, but af- 
ter inflation this term is irrelevant, and the second term 
in ([T]) is dominate after inflation. Therefore the energy 
stored in fields are to be transferred to other fields by 
thermal bath of fields. 

But if we demand that the energy density at the be- 
ginning of radiation epoch is the same as at the end of 
inflation, the thermal bath is not sufficient, due to the 
expansion of universe. So, a coupling between fields and 
matter is needed. 
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We suppose that the fields decay into relativistic parti- 
cles, X) with decay rate, T, which depends on details of 
interactions between the fields with the relativistic parti- 
cles. Here, we will obtain a bound on T from conservation 
of energy. We have [l[ 



Px + 3H (Px+Px) = T Pi>±- 



(35) 



Here 

Px/Px = V 3 ' Pl>-/Pl>- = W eff 



-1+gcW. (36) 



where w e ff is the equation of state. Hence, we can solve 
([3"5]) as 



Px 



(t) ujH' 2 (t r ) 7 wtr uj (J 



4(l)H*(l) e £(*v-'0 d 7 



(37) 

where t r is the time just after the end of inflation and M 
is the scale energy of H (t) that we explicitly show. After 
short time H — > l/2t, and a{t) — > dot*. By assuming 
that r»w, the fields almost immediately decay into x, 
i.e. H{t r ) w M, so 



~ A/>±(*r-) 



a(t) 



(38) 



Therefore, to have successful inflation and reheating with 
this scenario, we need HAti n f >• 1 and TAti n f ^> 1, one 
can set r ss H . For the standard scalar field, to produce 
a successful radiation epoch after reheating period, we 
must take r > H Q. 



VII. 



SUMMARY 



We have studied the Gauge-nation by analytic meth- 
ods and we have investigated the simple (but nontrivial) 
ansatz that shows the main features of the model. Then, 
we have derived formulas for leading order of fields in the 
model. The formulas are valid in all range of history of 
the early universe. 

Using the formulas, we have provided analytic solutions 
for the MS solutions [H, and with the analytic solutions, 
we studied some features of the MS solutions which can- 
not be obtained without analytic methods. Then, we 
obtained constraints from slow roll conditions on the pa- 
rameters of the solutions. 

Moreover, we studied preheating period in the model and 
obtained a bound on decay rate of fields, that may be 
useful for future works. 



ACKNOWLEDGMENTS 

I am grateful for helpful discussions with F. Arash, H. 
Asgari, M. M. Sheikh- Jabbari and A. Malcknejad. 



Appendix A 

In this letter we give some solutions for Eq. ([6]), It 
is necessary to show that they are also solutions of the 
equation of motion for fields. To show this point, in this 
appendix we will obtain the equation of motion for fields 
in terms of variables that we use in this letter. 
The equation of motion can be obtained by variation of 
(P) with respect to the fields as [1] 

{1 + K9 ^ } a +(1 + K ^ )29 ^ + ( 1 - 3 ^ ^ H a =°- 

(Al) 

But, another standard way, that we use here, to obtain 
(|A1[) is to use the Friedman equations. For what we will 
do, let us review this method. 
From Q, we obtain the following equations 




(A2) 



and 



(A3) 



the derivative of (|A2[) with respect to time results in 

„2 1,4,12 \ 



2HH = — - -Tr + ^—7- + n 2 „ 



2 dt 



(A4) 



H 2 = ^[ 



Substituting Eq. (|A3jl into the left hand side of (|A4[) , 
with algebraic manipulations, gives us (|AI[) . 
One way to obtain the equation of motion in terms of 
variables that we use in this letter, is to substitute vari- 
ables in Eq. (|AI|) . but it is better to rewrite Eqs. (|A2[) 
and (|A3[) in terms of our variables, and then to derive 
the equation of motion. 
For Eq. (|A"2j) . we have 

^-6 2 (t)H 2 sm 2 2a(t^j , (A5) 

and for Eq. (| A3|) . we have 

H=-e(t)H 2 . (A6) 

The derivative of (|A5[) with respect to time, substituting 
Eq. (|A6[) into the left hand side of the result, is 



-2H 3 e(t) = HH ( e[t) + ^-( 2 {t)H 2 sin 2 2a{t) 



+ \H 2 Ue(t) + ^e 2 (t)H 2 ^ 2 Mt) 



We rearrange Eq. (|A7J) as 



= -e{t)H 3 1-2 + e(t) + —e 2 (t)H 2 sin 2 2a(t) 
+ \H 2 j t Ut) + ^e 2 (i)H 2 s\n 2 2a(t) 



(A7) 



(A8) 
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Eq. (|A8[) is the equation of motion of fields in terms of then, from Eqs. (|A8|) and (|A9|) . we have 
our variables. 

If we rewrite Eq. © as _ e{t)H Jgg) + \± (L.H.S.fg®) = °- ( Al °) 



2 Therefore, solutions of Eq. (|A9[) are the solutions of Eq. 

e(i) + C*)^ 2 sin2 2a(*) - 2 = 0, (A9) (lAlO)) 
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